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A novel superconducting state under the broken time- reversal symmetry is studied in conven- 
tional phonon-mediated superconductors. By solving the Eliashberg equation self-consistently 
with the mass renormalization effect, it is found that the even- and odd-frequency components 
of the order parameter coexist in the bulk system as a consequence of the broken time-reversal 
symmetry. This finding would direct more attention to the odd-frequency pairing that affects 
physical quantities, especially in strong coupling superconductors. 
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Since superconductivity was discovered by Kamerlingh 
Onnes about a century ago, it has been one of the most 
important subjects in the condensed matter physics. 
Odd-frequency superconductivity is one of the issues at- 
tracting interest in this field today. It was originally sug- 
gested by Berezinskii in the "^Hc superfluid phase^ and 
was developed further by Kirkpatrick and Belitz^''^ and 
by Balatsky and co-workers^"^ in conduction electron 
systems. The possibility of the odd-frequency supercon- 
ductivity was also discussed in doped triangular anti- 
ferromagnets^ and in electron systems with critical spin 
fluctuations in the vicinity of the quantum critical point. ^ 
In a similar context, quasi-one-dimensional organic sys- 
tems were also studied.^' ^° 

In these works, the random phase approximation was 
frequently applied to examine the stability of the odd- 
frequency superconductivity, where the mass renormal- 
ization effect was neglected completely. However, it was 
pointed out that mass renormalization tends to desta- 
bilize the odd-frequency superconductivity.^ Thus, a 
more elaborate discussion on the realization of the odd- 
frequency superconductivity with the appropriate mass 
renormalization and vertex corrections is required. 

In this respect, the odd-frequency superconductivity 
is not so easy to emerge by itself. Alternatively, it is 
much easier to emerge or to be induced with the help of 
even-frequency superconductivity. This is analogous to 
the parity mixing of Cooper pairs, such as s-j-p-wavc, un- 
der the broken inversion symmetry. ^^ Similarly, the even- 
and odd-frequency components of the order parameter 
should coexist under the broken time-reversal symme- 
try. The possibility of such an induced odd-frequency su- 
perconductivity has also been examined in symmetry re- 
duction at interfaces^^"^^ or vortices, ^^ and in fully spin- 
polarized electron systems with orbital fluctuations.^^ 
Since the previous studies did not examine the self- 
consistent solution of the order parameters in the super- 
conducting phase, there is still room to reexamine the 
induced odd-frequency superconductivity. 



In this Letter, we explore the coexistence of the even- 
and odd-frequency pairings in conventional phonon- 
mediated bulk superconductors on the basis of the 
Eliashberg theory, where mass renormalization is appro- 
priately taken into account. ^^ We also emphasize the im- 
portance of an explicit self-consistent procedure of the 
Eliashberg equation in the presence of the odd-frequency 
pairing. 

As the simplest source breaking the time-reversal sym- 
metry, a uniform external magnetic field coupled to con- 
duction electron spins is introduced here. We note that 
other magnetic sources, e.g., ferromagnets or ferromag- 
netically polarized impurities, give essentially the same 
consequences. To stabilize superconductivity, we consider 
the electron-phonon interaction as usual. To make the 
discussion clearer, we focus on an Einstein phonon char- 
acterized by a single excitation energy. The electron- 
phonon interaction gives rise to attractive interactions 
for both spin-singlet and spin-triplet pairing channels. ^° 
The local nature of the Einstein phonon favors the s- 
wave pairing. Since the triplet s-wave pairing is forbid- 
den with the even-frequency dependence owing to the 
fermion property, we must extend the theoretical frame- 
work including the odd-frequency dependence.^' ^° 

We treat the electron-phonon interaction as the per- 
turbation and derive the lowest self-energy in 4 x 4 matrix 
form as follows: ^^'^^ 

S(itj/) = -g^-^T^D{iuji - i^„,) y^ p;jG{iuj,n, k)p^. 

m k 

(1) 



Here, g represents the coupling constant for the electron- 
phonon interaction. Np is the number of phonon sites. 
n represents the volume of the system. D([uji — iujm) 
and G(iLOm, k) are Green's functions for the phonon and 
conduction electron: 

Dim) ^ ^ 



wi - We 



-wi - We 



Giiui, k) = [iw; - efeP3 - hp^az - S(iw/)]" 



(2) 
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Here, we represents the phonon frequency, efc is the ki- 
netic energy of the conduction electron measured from 
the Fermi energy, h represents the external magnetic field 
along the z-axis. lyi = 2ttTI and cui = TrT{2l + 1) arc Mat- 
subara frequencies at the temperature T for bosons and 
fermions, respectively, p,^ and cTq, (a = 1,2,3) are Pauli 
matrices for the particle-hole and spin spaces, respec- 
tively. The self-energy S(iaj;) in cq. (2) is given by 



S(iw,) = -iS^ + Efpgcrg + EA(iwO- 



(3) 



The first and second terms arc for the normal parts that 
renormalize the Matsubara frequency and magnetic field, 
respectively. There is no p^ term in cq. (3), since it dis- 
appears after the summation over k in cq. (1) under the 
particle-hole symmetry. The third term in eq. (3) is the 
anomalous part. 

To consider the anomalous part, let us discuss what 
types of matrices come out from G{iuji,k) in eq. (2). 
Since the singlet pairing (p2'''2)^^ is stabilized at low 
temperatures, we assume that it is finite. Then, the prod- 
uct of the singlet component {p20'2) a-nd the Zeeman 
splitting term {p^cr^) gives rise to the following addi- 
tional matrix component: p-^(T:iP2cr2 oc Px""!- Moreover, 
a finite Picri component results in p^crspicri oc P2T2. 
The PiCTi component represents the Sz ~ spin-triplet 
pairing. Since the orbital of the Cooper pair is of the 
s-wave type, the Picri component must exhibit the odd- 
frequency dependence.^ Therefore, the anomalous self- 
energy should have the following form: 



SA(iw;) = Afp20-2 + iA°PiCri. 



(4) 



Here, the first term A°p2cr2 represents the s-wavc singlet 
pairing with the even-frequency dependence. The second 
term iA°picri is for the s-wave triplet pairing with the 
odd-frequency dependence. With this choice of the rela- 
tive phase between the even- and odd-frequency compo- 
nents, A^ and A° can be chosen as real quantities in the 
Eliashberg equation, as shown later. 

In the presence of the magnetic field, the Sz value along 
the field is a good quantum number. However, the total 
S is not the case. Therefore, the singlet and Sz ~ 
triplet components are mixed, while the other Sz = ±1 
triplet components are decoupled completely from the 
singlet component. ^^ Thus, the s-wave triplet with the 
odd-frequency dependence appears inevitably with the 
singlet under the magnetic field. In fact, this type of self- 
energy was taken into account years ago to study the pair 
breaking effect caused by magnetic fields or magnetic 
impurities. ^'^ At that time, however, the self-consistent 
treatment for A° was not performed, since it was prema- 
ture to consider the odd-frequency superconductivity. 

In this work, we solve the Eliashberg equation self- 
consistently and reveal that the odd-frequency compo- 
nent appears simultaneously with the even-frequency one 
below the superconducting transition temperature Tc. 
Substituting eq. (2) into eq. (1), we obtain the Eliashberg 
equation in the following matrix form: 






^xtj: 



m>0 



vr f+ 
V,'^ f- 



-V,- 



' ImJm 
Im-J m 



UJr, 



A? 



A^E : 



r?i>0 



V+ f+ 

-vr f~ 

ImJ m 



Im-J m 
V,~ f+ 



A^ 
A° 



(5b) 

Here, A ~ 2g^A^pA^(0)/ajEO is a dimcnsionless coupling 
constant with A^(0) as the density of states per volume 
at the Fermi energy. w,„ and ft,,„, which are the renor- 
malized Matsubara frequency and magnetic field, respec- 
tively, are as follows: 

Lbyn=UJ.m+T!:^^, hni^h + Y.^^. (6) 

^irn ^'"^ effective interactions for the even-frequency (+) 
and odd- frequency (— ) components with respect to both 
the LJi and ujm frequencies. They are given by 



V, 



Irn 



± 



(7) 



[ioi - a;,„)2 + w| (cj; + w„i)2 + a;! 

/j^ used in eqs. (5a) and (5b) are defined as the real and 
imaginary parts of 

■K 



^™ + i/i™ +(-A;;, + iA^)' 



(5a) 



^I7a + ^f;n- (8) 



Here, we adopted a sufficiently large cutoff in the Mat- 
subara summation and checked that there is no cutoff 
dependence in the self-consistent solutions. 

The two matrix equations in eq. (5) are coupled via 
the fi^ functions that contain the normal and anoma- 
lous parts of self-energies. We note that the superscripts 
± in f^ represent even and odd functions as in Vj^. 
Therefore, eq. (5) clearly shows that Sj* and A^ (EJ^ and 
A°) have the even (odd) -frequency dependence. We note 
that the relative sign between A° and A® is the same as 
the sign of /i. In the absence of the external field {h = 0), 
it is expected that Am = 0, /~ = 0, and AJ^ = 0. In 
this case, eq. (5) reduces to its typical form for the even- 
frequency superconductivity. 

First, we determine Tc by solving the linearized gap 
equation. For this purpose, we put A^ = A°„ = in cq. 
(8) and calculate the normal self-energy using eq. (5a). 
With the use of the determined Cum and /i,„ , the gap equa- 
tion eq. (5b) reduces to an eigenvalue problem. Tc is de- 
termined when the eigenvalue is equal to 1. For /i — 0, the 
matrix for the eigenvalue problem reduces to ATVJ^/+. 
In this case, all the eigenvalues are real, since the matrix 
can be symmetrized. In contrast to this, the matrix is 
no longer Hermite for /i 7^ owing to the odd-frequency 
component and the eigenvalues are complex in general. In 
Fig. 1(a), we show the largest two real eigenvalues under 
various fields, where the coupling constant is chosen to be 
as relatively strong as A = 10 to emphasize the effects of 
the odd-frequency component. For h ^ 0, the real eigen- 
values exist only in the low-temperature region. We can 
see that Tc decreases with h. In Fig. 1(b), we show the 
phase boundary on the T-h plane. To see how the odd- 
frequency component affects the stability, we also exam- 
ine Tc for a pure even-frequency case by solving eqs. (5a) 
and (5b) self-consistcntly with A°^ = 0. This is exactly 
the solution obtained on the basis of the traditional the- 
ory only within the even- frequency component. We note 
that Tc is suppressed by the induced odd-frequency com- 



J. Phys. Soc. Jpn. 



Letter 



M. Matsumoto, M. Koga, and H. Kusunose 



2 
1.8 
1.6 
1.4 
1.2 

1 
0.8 
0.6 



\ (a) : 




h=1.5 

h=2.3 
h=2.8 


N. \ : 


\ \ : 


h=3.0 " 


_ :>», ^>s,.;_^^ 


-— __i.______- 




1 


1 1 1 



100 



0.1 



0.2 



0.3 



0.4 

T/(Bp 



0.5 



0.6 



0.7 




T/(Op 



0.1 



! 

(b) ____:___ 


i even+odd 

even 




i '"^v^ 




^^^' '■' 





0.2 



0.3 



0.4 



0.5 



0.6 



12 
10 

. 8 

3 

E 6 

3 

' 4 
2 






\ X ^ h=2.3 (even) 

- \ X : h=2.3 (odd) 

\ \ h=2.8 (even) 

\ \ h=2.8 (odd) 


\ : ^\ : T=0.1o)e 


---K.^\ 1 


1 1 1 ~- — 1 ----- 



10 



T/(Of 



CB^/COp 



Fig. 1. (Color online) (a) T dependence of the largest two real 
eigenvalues for A = 10. The magnetic fields h are given in units 
of cjg. (b) Phase boundary on the h-T plane. The solid line 
represents the results of the coupled even- and odd-frequency 
components, while the dashed line is obtained by retaining the 
even-frequency component alone. 



poncnt. Since Tc(cvcn) is higher than Tc (even + odd), 
one may think that the pure even-frequency pairing is 
more stable and that the odd-frequency component does 
not appear. However, the Ehashberg equation (5) clearly 
indicates that the even-frequency component alone can- 
not be the solution under a finite field. We emphasize 
that the appearance of the odd-frequency component is 
inevitable in the absence of the time-reversal symmetry 
as h j^O. 

Next, we examine the self-consistent solution of the 
Eliashberg equation. Figure 2(a) shows the temperature 
dependence of the order parameters that arc summed 
over the positive Matsubara frequencies to indicate the 
strength of the frequency-dependent order parameters. 
We can sec that the odd-frequency component appears 
below Tc for /i 7^ and that both components develop at 
low temperatures. The odd-frequency component is en- 
hanced by h and the two components become comparable 
under high fields, since the off-diagonal matrix element 
ifm) becomes dominant in eq. (5). In an extremely high 
field case, however, inhomogeneous solutions, such as the 
FFLO state solution, should be taken into account. ^^'^^ 
In contrast to the FFLO state, we emphasize that the 
effect of the odd-frequency component is present even 
under low fields. Figure 2(b) shows the frequency depen- 
dence of the order parameters for a fixed temperature. 
The even-frequency component has a maximum value at 



Fig. 2. (Color online) (a) T dependence of the even- and odd- 
frequency components under various fields for A = 10. The 
strength is defined by 5Zm>0 ^(''^'")/'^E- (b) i^m dependence 
of the order parameters at T = O.Ioje- 



the lowest Matsubara frequency, while the odd-frequency 
component forms a peak structure at cjm — ^e- This is 
because the attraction for the odd-frequency component 
vanishes at low energies, as indicated by eq. (7). 

For a weaker coupling constant A, Tc shifts toward the 
lower temperature region as usual. The induced odd- 
frequency component decreases accordingly. Even in a 
weak coupling case, however, the odd-frequency compo- 
nent appears below Tc as long as the magnetic field is 
finite. As in Fig. 2(a), it is enhanced with the increase in 
magnetic field. 

Finally, we elucidate the stability problem of the gap 
equation, which is closely related to the long-standing 
controversy in the odd-frequency pairing. That is, it was 
argued that the pure odd-frequency superconductivity 
was accompanied by a negative Meissner kernel^ and 
was thermodynamically unstable in the bulk.^^ How- 
ever, the path-integral formulation^^ revealed that the 
odd-frequency superconductivity shows a typical Meiss- 
ner effect when the effective free energy is minimized 
appropriately.^®' ^^ We discuss this problem in view of 
gap equation solvability. For this purpose, we rewrite the 
anomalous self-energy in a generic form by introducing 
the explicit global phase factor e"''. In the coexisting case, 
the gap equation requires the following form: 



SA(iw/) 




,(P) 
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,(P) 



.(P) 



Aze" 



.(h) 



.(h) 



-A? 



(9) 



iA° 



iA° 



Here, AJ and AJ are the order parameters for the 
particle and hole components, respectively. The relative 
phase is fixed as " i " with the real A^ and A° as it should 
be from the structure of the gap equation. To maintain 
the common global phase, the factor (f> must be chosen 
as either (/) = +! or 0=—l. Then, the denominator in 
eq. (8) is altered accordingly as 



\l iibrn + i/lm j 



'(-A^ + iA5^)' 



(10) 



It should be emphasized that the solution of the gap 
equation can be found only when we choose the correct 
sign of (j). Namely, from the structure of the gap equation 
[eq. (5)], the following condition in eq. (10) is necessary 
to obtain the nonzero solution: 



'(-A°„ 



iA°j'>0. 



(11) 



This is an important clue to understanding the stability 
problem. 

Along this line, let us discuss first the case of a pure 
even- or odd-frequency superconductivity. When h = 0, 
a pure even-frequency superconductivity is realized. In 
this case, the solution of the gap equation can be ob- 
tained only when we choose ^ = -|-1 in eq. (10) as usual. 
On the other hand, in the case of a pure odd-frequency 
superconductivity, which can be realized by neglecting 
the mass renormalization in eq. (5),^*^ the solution can 
be found only when we choose (j) = ~1 owing to the 
prefactor " i " of A°„ in eq. (10).^° When we choose the 
wrong sign, (j) = +1, no solutions can be obtained for the 
pure odd-frequency superconductivity.'^'' In summary for 
the pure even- or odd-frequency superconductivity, the 
solvability condition of the gap equation is represented 

by 



.(p) 



t 



.(h) 



(12) 



This condition is consistent with that proposed by the 
path-integral formulation and the minimization of the 
free energy.^^"^^ 

Next, we discuss the case in which the even- and odd- 
frequency superconductivities coexist. In this case, Af 
(A°) is the majority (minority) component. Here, A° is 
induced by the magnetic field maintaining the relative 
phase " i ". According to the key condition, eq. (11), the 
majority component between A^ and A° determines the 
correct sign of (j). That is, the even-frequency component 
is the majority component in the present problem and 
gives the sign (p — +1. Indeed, the solution of the gap 
equation can be found only with ~ +1. In the coex- 
istence case, eq. (12) is no longer satisfied. It is valid 
only for the pure even- or odd-frequency superconduc- 
tivity.^^"^^ Instead, we can generalize the stability and 
the solvability condition for cj) so as to satisfy eq. (11). In 
any cases, we note that the generalized condition (11) is 
consistent with the minimization of the free energy. This 
will be discussed elsewhere. 



In summary, we studied the superconducting state 
under the broken time-reversal symmetry in conven- 
tional phonon-mediated superconductors. We solved the 
Eliashberg equation self-consistently with a proper mass 
renormalization. It is found that the even- and odd- 
frequency components of the order parameter coexist as 
the consequence of the broken time-reversal symmetry. 
The induced odd-frequency component affects more or 
less various physical quantities, such as transition tem- 
perature, density of states, Mcissncr effect, and Knight 
shift. In the Meissner effect, for instance, the induced 
odd-frequency component gives a paramagnetic contri- 
bution to suppress the Meissner screening. This is be- 
cause the minority component suppresses the majority 
component owing to the relative phase factor " i " be- 
tween A°„ and A°„ in eq. (8). It should be noted that 
the paramagnetic Meissner effect never exceeds the typ- 
ical diamagnetic Meissner effect, since the former comes 
from the minority component. The role of such an odd- 
frequency component of the order parameter has not 
been considered in past studies on superconductivity. It 
is likely prominent in strong coupling superconductors, 
such as Pb, where a considerable enhancement of the in- 
duced odd-frequency component could be expected. It is 
worth seeking for such effects experimentally under mag- 
netic fields or in ferromagnetic superconductors. 
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